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Abstract 

In this paper we calculate the pressure of pure lattice Yang- Mills theories and lattice QCD with 
heavy quarks by means of strong coupling expansions. Dynamical fermions are introduced with 
a hopping parameter expansion, which also allows for the incorporation of finite quark chemical 
potential. We show that in leading orders the results are in full agreement with expectations from 
the hadron resonance gas model, thus validating it with a first principles calculation. For pure 
Yang-Mills theories we obtain the corresponding ideal glueball gas, in QCD with heavy quarks our 
result equals that of an ideal gas of mesons and baryons. Another finding is that the Yang-Mills 
pressure in the large N limit is of order ~ N° to the calculated orders, when the inverse 't Hooft 
coupling is used as expansion parameter. This property is expected in the confined phase, where 
our calculations take place. 



1 Introduction 



With the experimental heavy ion programmes at RHIC and LHC, it is possible to perform experi- 
mental research on the properties of the quark gluon plasma. In order to distinguish this new matter 
phase from ordinary matter, it is decisive to have an accurate knowledge of both. While in the decon- 
fined phase qualitative insights are often based on perturbation theory, one of the main theoretical 
frameworks used to describe the hadronic phase, besides e.g. chiral perturbation theory, is the hadron 
resonance gas model and its variants. The idea is to model the confined phase as an ideal gas of 
hadrons and resonances. Confronting predictions based on this picture with non-perturbative results 
from lattice simulations shows good agreement up to the phase transition region [HI]- 

At present, the hadron resonance gas is a model requiring a lot of experimental input data on 
hadronic properties in order to give predictions, and it lacks a proper field theoretic motivation. In 
this paper we derive the hadron resonance gas model as an appropriate effective theory for strong 
gauge couplings from the QCD Lagrangian with quarks and gluons as the fundamental degrees of 
freedom. In particular, we calculate the pressure - or the negative free energy density - in a combined 
strong coupling and hopping parameter expansion in the lattice theory. These expansions are known 
to give convergent series. Strong coupling methods have also been used, e.g., to employ the properties 
of thermal QCD with mean field methods [MI], numerical simulations [8] or series expansion methods 

[SHI]. 

The outline of this paper is the following: after a short review of the hadron resonance gas model 
and its description on the lattice, we calculate the pressure for pure SU(N) lattice Yang-Mills theories, 
with special emphasis on the cases N = 3 and N — > oo. Then we introduce dynamical quarks via a 
hopping parameter expansion and calculate the pressure for QCD with heavy quarks. In all cases we 
obtain the right spin degeneracies for the lowest mass hadrons and the correct dependence on chemical 
potential when compared to predictions of the hadron resonance gas model. 



2 Hadron resonance gas 

The basic quantity of thermal systems is the grand partition function [12j 

Z = tre-<*-" JV > T "\ (1) 
from which we can deduce the pressure as 

dV v^oc V y ' 

The last equality holds in the infinite volume limit, which we will always consider in the following. 

In the ideal gas picture, we can write the partition function as a sum over one-particle partition 
functions of the hadrons and resonances [T][5] 

\nZ = J2^Zl (3) 

i 

In the continuum, the one-particle partition functions at temperature T and chemical potential fi{ are 
given by 



oo oo oo 



\nZl = V 9i J J /^ln(l + ,e— ), ( 4 ) 



-oo — oo — oo 
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where gi is a possible spin-isospin degneracy factor and w$ = \/p 2 + mf for particles of mass rrii. The 
factor 77 refers to bosons (77 = —1) or fermions (77 = 1). With Eqs @ and ([3]) the pressure of the 
hadron resonance gas is given as a sum over one-particle contributions 



P 



(5) 



Calculating the ideal gas pressure on the lattice, we have to keep in mind that the momentum inte- 
gration is limited to the first Brillouin zone, i.e. to the interval [ — — , — ). Hence the pressure is given 
by 



Pi = 



d 3 P 



77 In 



1 + 77c 



-{ui—iii)L T 



where the temperature is related to the temporal lattice extent via 

T = 1 1 



(G) 



(7) 



For heavy particles it is convenient to use non-relativistic approximations. In the static limit we have 
to — ttl. However, in the continuum this approximation leads to a divergent momentum integration. 

Moreover, we have to integrate over arbitrarily large momenta where the approximation does not hold. 

2 

Thus we have to expand the relativistic dispersion relation at least to ui ~ m + f— in order to render 
the integral finite. 

Using the static approximation on the lattice, the momentum integration gets trivial and we obtain 



Pi 



9i 



-77 In 



f 



?7e 



-(mi— fii)L T 



(8) 



In the following we will use the hopping parameter expansion, which is applicable for heavy quarks 
and thus for heavy composite particles like the corresponding mesons and baryons. Hence it is justified 
to consider only the limit 



L T 



which yields 



P 



P, 



-(nii-iii)L T 



' L T a 3 

Using the fact that to each baryon with fis there is an antibaryon with —fJ-B, we obtain 
P = t~ — t V g M e- mML - + V g B e- mBL ^ cosh(/i B L T ), 



L T a 3 



(9) 



(10) 



(11) 



B 



where we have split the sum into a mesonic and a baryonic part. It is this form of the pressure which 
we will derive in the following from a first principles strong coupling calculation. 



3 Pressure of lattice Yang-Mills theories 

In an earlier paper [11] the general procedure to calculate the pressure in case of SU (2) has already 
been discussed. Thus we will only shortly outline the derivation of the SU (N > 3) results and focus 
on the more interesting cases N — 3 and the limit N — > 00. 
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Figure 1: Graphs that appear on the different iattices. Left: Leading order tube for N T — 4. Middle: 
Graph with an inner plaquette in the case N T — 4. Right: Graph with two slits contributing on both, 
N T = 4, oo, lattices. 



To obtain the physical pressure, we need to calculate 

Pphy S = P(ir)-PM, (12) 

where the subtraction of vacuum contributions (L T — > oo) is necessary for renormalisation. On a 
lattice with temporal extent L T the expression to evaluate is 

P(L T ) = ^lnc (/3) + -3- a^Il^r- ( 13 ) 

C=(X" 1 ) * 

Here &(Xi) is the contribution of a graph Xi, C specifies a cluster of graphs and a(C) is a combinatorial 
factor determined by the formalism of moments and cumulants (13) . When L T goes to infinity this 
is the zero temperature contribution. The leading order contribution to the subtracted pressure is 
a tube of AN r fundamental plaquettes spanning around the temporal direction, as shown in Fig. [T] 
(left). In contrast to SU(2) we have two fundamental representations in the character expansion and 
counting the 3 different cross sections of the tube per site we get 

where u is the character of the fundamental representation and effective expansion parameter. A 
volume factor cancels out, due to the spatial translations of the tube. This is the leading order 
result for all N > 3. An explicit dependence on the gauge group appears in higher orders, but is 
indirectly encoded in the functional behaviour 



»(■>'> = 2^a + ---- lVt] 



It is interesting to note that for N — > oo we have the correspondence [2] 

^)-^ = ^=A- (16) 

and hence the expansion parameter of the fundamental representation equals the inverse of the 't 
Hooft coupling A [15] . 

The leading corrections to the basic polymer arise from putting spacelike plaquettes in higher 
dimensional representations inside the torus, as in Fig. [1] (middle). In order to have non- vanishing 
contribution according to the integration formula 

dU Xr {U) = 5 rfi , (17) 

the allowed representations are those whose conjugates appear either in / ®/ or in /®/, the Kronecker 
products of the fundamental representations. The latter representations are forced to appear an even 
number of times along the tube, since they induce a change of orientation of the outer plaquettes and 
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Figure 2: Orientation of plaquettes that meet at some link (thick line) in case of SU(3). Left: 
Fundamental and antifundamental representation at the link resulting in a conserved orientation of 
the corresponding plaquettes. Right: Three fundamental representations at the link and a change of 
orientation, which has to be cancelled to give an allowed graph. 



this has to be cancelled due to the periodic boundary conditions, cf. Fig. [5J If orientation changing 
plaquettes appear n times, we take this fact into account by inserting a factor 

\ [! + (-!)"] (18) 

into the summation of these contributions. Let us denote the parameters of the representations 
appearing in / (g) / with v\ and v 2 and their respective dimensions d Vl , d V2 . In general, 

d VlVl = N{N 2 + l) (Nuf + ... (19) 

d V2 v 2 = N[N ~ l) v 2 = \ (Nuf + .... (20) 

In the case of 5[/(3) we have to identify v 2 with the complex conjugate fundamental representation 
(3 €5 3 = 3 © 6), hence in this case v 2 equals u rather than being of 0(u 2 ). The parameter of the 
adjoint representation appearing in / ® /, which we will call w, fulfills 

d ad w = (N 2 - 1) w = (Nuf + . . . . (21) 

Another correction is to add slits to the tube by inserting two fundamental plaquettes at the same 
location, cf. Fig. [T] (right). For finite L T the minimal number of slits is two, since no plaquette is 
allowed to be occupied twice by the same polymer. This insertion gives an additional factor of 27V 2 u 2 . 
The factor of two appears because adding a slit increases the number of polymers by one and this 
new polymer can contribute in two different orientations. Thus we get a chain of polymers with a 
combinatorial factor a(C) of the resulting cluster C, which depends on whether the chain polymer is 
embedded in the finite lattice with periodic boundary conditions or in the infinite lattice. It is given 

by 

Finite lattice, pbc: a(C) = (i - l)(-l) i+1 

Infinite lattice: a(C) = (-if, (22) 

where i counts the number of polymers. Now we sum over all possible contributions of the torus with 
inner decorations. This calculation proceeds in the same way as in the case of SU(2) [llj . with the 
difference of inserting the factor Eq. (|T8|) . The result splits into two terms, 



P = 



—u 4N - 



c Nt + b N ^ 



(l + 0( u 4 )). (23) 



We have introduced 



cat = l + d Vl vi+d V2 v 2 + d ad w-2d 2 u 2 (24) 
b N = 1 - d Vl vi - d V2 v 2 + d ad w (25) 
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Figure 3: Correction graphs. Left: One shifted plaquette. Middle: One added cube. Right: Graph 
with a perimeter of six. To all graphs one can also add inner plaquettes and slits. 



as gauge group dependent factors. In the limit N — »■ oo the expansion parameters simplify [2] 



d Vl vi = -(Nu) 2 , 



d V2 v 2 = 

dadW = 



(Nu) 2 , 



1 

r 

(Nu) 2 , 



(26) 



and these gauge group factors become trivial 

1 



2„,2 



1 + -(Nu) 2 + -{Nil) 2 + (Nil) 2 - 2N 2 u 



6oo = 1 - -(Nu) 2 - -(Nu) 2 + (NuY 



1. 



(27) 



As a consequence the contributions of inner plaquettes cancel against each other. This gives a first hint 
that in the confined phase the pressure is a quantity of O(N ), and thus finite in the limit N — > oo, 
when expressed in terms of the inverse 't Hooft coupling u. 

It is in complete analogy to ST/ (2) to calculate higher order corrections. One can decorate the 
leading polymer, either geometrically or using higher dimensional plaquettes, one can add whole new 
polymers or one can combine these possibilities. Examples are shown in Fig. [3l Again we have 
contributions of some polymers only for large N T , e.g. the graph in Fig. [3] (left), which is absent for 
N T = 1. Some polymers only contribute for small N T , since their correction is of order u 2Nt or higher, 
like the graph in Fig. [3] (right). The results for the series up to 0(u a ) in the correction for the cases 
N = 3 and N — > oo are given in the appendix. 

3.1 Hadron resonance gas for SU(3) 

In the case of SU(3) the three lowest lying glueball states, which can be extracted from the strong 
coupling series of plaquette correlators, have masses of [TTHTg] 

ami (A) = -4 In u - In (l + 3u + 6vi + 8w - 18u 2 ) + 0(u 4 ) (28) 
am 2 (E) = -41nu-ln(l + 3w + 6vi+8w-18u 2 ) + Q(u A ) (29) 
am 3 (T) = -4 In u - ln(l - 3u - 6v 1 + 8w) + 0(u A ). (30) 

The subscript denotes the spin degeneracy of the glueball. Rewriting our result of the pressure 

3 



P 



4N T 



L T a 3 



■Mr 



(l + 0(u*)) 



we obtain 



L T a 3 



c -yn 1 (A)L T . 2 e ~ m 2(-E)L T _j_ g e -m 3 (T)L T 



(l + 0(u*)) 



(31) 



(32) 



which has to be compared with Eq. (fTTj) . We see that the hadron resonance gas result equals the one 
derived with strong coupling methods to the orders considered here. Corrections to the ideal hadron 
gas results appear in 0(u 4 ). 
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Figure 4: Polymers contributing to O (u 6 ) in the correction. Left: One cube has been added inside the 
tube. Middle: Four fundamental plaquettes have been changed for plaquettes ~ u 2 . Two fundamental 
plaquettes had to be added at the inside. Right: A graph with temporal extent N T — 5, not fitting 
into the compactihed lattice and thus contributing with a negative sign on the infinite lattice. 



The strong coupling expansion of the pressure can also be used for the determination of the lower 
integration constant in the so-called integral method. In numerical simulations a direct computation 
of the partition function is not possible jTB] , so that one has to calculate derivatives to get physical 
quantities. The pressure can be calculated, e.g. as an /3-integral 







* f dp'[6-p -3(V T + V s )], (33) 



where V s ,t are the expectation values of space-space and space-time plaquettes on an asymmetric lat- 
tice N T < N s and Vo is the one of a symmetric lattice N T = N s . The lower integration constant P(Po) 
has to be fixed by hand. Usually one argues with an exponentially small pressure P ~ exp(— ma/T) 
below T c , where ma is the lowest glueball mass, and chooses /?o correspondingly. With our calcula- 
tions we have justified this assumption, but it is also possible to use our series expansion to consider 
some specific value P(/3q) for larger values of /3o- 

3.2 The limit N -> oo 



We have already seen that to leading order in the strong coupling expansion the pressure is of O(N ) 
in the confined phase. Now we will give some examples of higher orders, where cancellation of all 
iV-dependence takes place. Let us consider at first the leading geometric decoration, a shift of one 
plaquette of the tube, resulting in an 0(u AN ^ +A ) graph, like the one in Fig. [3] (left). In place of the 
shifted plaquette, we can add v± , v 2 and w a d plaquettes weighted with their respective dimensionality. 
We can also put a completely new cube at this place, coming with a factor —2N 2 u 6 . A graph of this 
type is shown in Fig. [3] (middle). All possibilities summed up, we obtain the contribution 



P = 



'n 



I + 12AT T (l + d Vl Vi + d V2 v 2 + d ad w - 2d 2 f u 2 )u 4 + O (u 6 ) 



(34) 



The factor of 12N T emerges because we can shift each of the 4A r plaquettes in three directions. In 
round brackets we recognize another term cn, which again cancels in the limit N —¥ oo, and we obtain 



P 



L T a 3 



1 + l2N T u 4 + O (u 6 



(35) 



where every iV-dependence has disappeared. 

Consider next the graphs of Fig. 01 defined on the N T lattice or the infinite lattice. They are of 



6 



order 0(u e ) in the correction to the basic polymer and their contributions are 



left: $1 
middle: $2 



L T a 3 
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right: $3 



L T a 3 
6 

L T a 3 
6 

L T a 3 
6 

L T a 3 



,4JV T 



,4JV T 



,4iV T 



2N T d 2 f u 6 



N T d)u~ 2 w A J , 



N T dju 6 



Noting that in the limit JV->oowe have 



we see that 



w — V\ — V2 = u , 



$ = $1 + $1 + $3 = 0. 



(36) 



(37) 



(38) 



The contributions of the left and the right graphs cancel those of the graph in the middle and the net 
contribution of these graphs vanishes. We remark that by leaving the limit, these graphs contribute 
with \/N corrections. 

These cancellations occur in similar ways up to O (u 8 ) to which we have calculated, but in higher 
orders they become more subtle. The overall picture we find is the following: fix some certain inner 
plaquette in a graph. The representation of this plaquette be of dimension ~ N q . Now draw all 
allowed graphs, possibly consisting of several polymers, where the product of all plaquettes at these 
original plaquettes places is also ~ N q . If we have calculated the correct contribution and taken the 
combinatorial factor a(C) properly into account, the TV-dependence cancels out. To conclude, judging 
from our truncated series expansions, we consider the pressure to be of O \N°) in the confined phase. 



4 Pressure of lattice QCD with heavy quarks 

After having calculated the leading orders of the pressure for pure lattice Yang-Mills theories, we 
want to introduce fermions in our calculations. The easiest way to do this is to employ the hopping 
parameter expansion. This is an expansion in k = (2am + 8) _1 and thus a good approximation for 
heavy quarks. 

4.1 Hopping parameter expansion 

As the fermionic action is bilinear in the quark fields it is possible to integrate them out exactly, which 
we will do for Wilson fermions. Performing the Grassmann integration |13j . we find for an arbitrary 
flavour 

-S* = -lndet(Q / ) = -trln (Q f ), (39) 
where the quark matrix is given by 

Q f = l-K f M, (40) 
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and the hopping matrix M by 




(41) 



This leads to the hopping parameter expansion by expanding the logarithm 



= -tr In (1 - Kf M) = X] ~T trM '" 



(42) 



The delta functions in the hopping matrix M force the sum to extend solely over closed loops on the 
lattice. 

For sufficiently strong coupling, we can omit the gluonic contribution to the partition function and 
write 



I f 1=1 

In principle the sum extends over all closed loops I on the lattice, which for finite lattice extents can 
have non-trivial winding numbers. In the finite temperature case, these are e.g. the Polyakov loops, 



Due to the antiperiodic boundary conditions of fcrmions, this gives an additional factor (—1)" to 
Eq. (|43p . which depends on the winding number n. 

One might wonder that for temporal lattice sizes N T > 4 the leading contributions to the effective 
action are those coming from the plaquette ~ k 4 . However, such contributions cancel in leading orders 
of the pressure when subtracting the zero temperature part. The leading terms are therefore those 
from the smallest graphs with a non-trivial winding number, the Polyakov loops. These terms are 
~ k Nt in the effective action, or ~ (ne ±afJ -) NT with finite chemical potential and the sign depending 
on the orientation of the Polyakov loop. In the confined phase the Polyakov loop cannot be screened, 
so the leading term in the pressure will consist of at least two (oppositely oriented) Polyakov loops, 
giving a factor of ~ k 2Nt . This corresponds to a mesonic contribution since we used two static quark 
lines, one for a quark and one for an antiquark. To get a full hadron resonance gas picture, we also 
need baryons, which consist of three quarks or antiquarks. Hence we also have to take into account 
generalized Polyakov loops that wind around the temporal direction two and three times before being 
traced. If we do so, the effective action leading to the lowest (physical) mass hadrons reads 




(43) 



L s = trW s = ti Y[U (x,t). 



(44) 




+ 




X 




(45) 



X 



The effective partition function is given by 




/ 



(46) 
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Expanding all terms up to 0(k 3Nt ) and doing the group integrals we get the following result for the 
pressure for two flavours up and down 



P = 




(47) 



If we now use the expressions for the hadron masses given in the appendix, we are able to rewrite this 

as 



where the sum over the pseudoscalar mesons includes the pions as well as the 77° and the vector mesons 
include also the cj°. This is due to the fact that we calculate with heavy quark masses and there is 
no reason for the pions to be lighter. We are far away from the chiral limit. The prefactors for the 
baryons include the spin degeneracy as well as a factor 2 counting the corresponding antibaryons. 
Quark chemical potential fi has been replaced by baryon chemical potential [is = 3/U. 

As in the case of pure Yang-Mills theories, we see that the results expected from the hadron 
resonance gas have been reproduced by our strong coupling and hopping parameter expansion. Similar 
formulas hold for a larger number of flavours. We hence conclude that at least in this parameter regime 
the hadron resonance gas model has been fully validated. 

5 Discussion 

In this note we have calculated the pressure of lattice Yang-Mills theories and QCD with heavy quarks 
in combined strong coupling and hopping parameter expansions in leading orders. Chemical potential 
has been introduced as an additional parameter. The leading order results are valid for all N T in the 
confined phase for small enough (3. They can be used to fix the lower integration constant P(f3 « 0) 
in the integral method. For larger (3 one should use the corresponding series expansion. In pure 
Yang-Mills theories we were able to calculate these expansions up to eight orders in the correction to 
the leading order term. In this case we could validate the hadron resonance gase model and indicate 
on the finitencss of the pressure in the large- N limit. 

In the case of QCD with heavy quarks we could also derive expressions that equal results of the 
hadron resonance gas model, albeit only in the strong coupling limit and leading orders in n. Future 
work could account for (3 and k corrections in this case. In this way one could study the effects of 
mass corrections, the dispersion relation as well as deviations of the ideal gas behaviour by hadron 
interactions. 
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P 




(48) 
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A Hadron masses from hopping parameter expansions 

In order to extract hadron masses, one has to specify suitable operators, which have non-vanishing 
overlap with the state, whose mass one wants to calculate. An interpolating operator for mesons A{x) 
can be constructed from a quark and an antiquark and some 7 matrices to specify the spin structure, 
like for the pions and p-mesons [13] 

7T + (x) = d ac (x) (l5) a pU0c(x), 

7r_(x) = u ac (x) (j 5 ) a pdp c (x), 
p+(x) = d ac (x)(j k ) a/3Ul3c (x), 

p-(x) = u ac (x) (7fe) ai g dp c (x). (49) 

The generalization to mesons with a different quark content and other vector mesons should be obvious. 

For baryons or antibaryons we can do the same using three quarks or antiquarks and some 7 matrix, 
e.g. proton and antiproton can be constructed via 

p%{x) = e cde {C-i 5 ) 01 u ac {x)[u l3d (x)d ie (x) -~ df3d{x)u ie (x)\, 
Psiv) = e fgh(C^5) ev u S f(y)[d Eg (y)u lph (y)-u Eg {y)d !ph (y)], (50) 
with the charge conjugation Dirac matrix C, satisfying 

C 7/i C- 1 = -7j, -C = C T = C'- 1 = Cl (51) 

Now we are able to calculate the masses of the lightest particles with these quantum numbers from 
the exponential decay of this interpolating operator in the following correlation function 

m(A) = — lim — In (A(0)AUt)) for mesons, 

t— >oo t 

m(B) = - lim -ln(B(Q)B(t)) for baryons. (52) 

t— >oo t 

In order to have mass eigenstates, one has to sum over all spacelike coordinates 

A(t) = ^A(S,t), (53) 

X 

but as we only want to calculate the leading order contribution, we need not care about this. 
The masses that we obtain in this way are quite simple. They are given as 

Mesons: am ff' = ~ \n2nf — \n2Kf> (54) 

Baryons: am//'/» = — \n2nf — ln2K// — ln2«/" (55) 
Antibaryons: amfppi = — \jx2kj — \n2Kp — In2/Cyr«. (56) 

We see that in leading order hopping expansion hadron masses are linear in the corresponding con- 
stituent quark masses with flavour /, if we define 

am/ = —]ja.2Kf O ft/ = — exp(— arrif), (57) 
relations that have been known for a long time [T9"ll20| . 
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B Strong coupling series of the Yang-Mills pressure 

The results for the pressure P(N T , u) of SU(3) up to higher orders and small N T are 

. \ 4 „ 4 c 7 2251587 o 641763 „ 

P(l,w)a 4 = 6u 4 + 44u 6 -30u 7 + u 8 u 9 

v ; 4096 1024 

124533967 10 
+ 20480 U ' 

n 4 o 8 in ii H05 12 351 n 18058839 14 

P(2,u)a 4 = 3u 8 + 27w 10 -135?i U + u 12 u 13 + u 14 

4 4 10240 

14775291 648558969807 1R 
2048 83886080 

t - * / o \ 4 „ i2 .,i4 is o , ifi 1269 17 199458821 1S 
P(3,u)a 4 = 2u 12 + 54u 14 -216u 15 + 315u 16 — u 17 + m 18 

89855001 19 671236701 20 
5120 U + 40960 U ' 

x 4 3 1fi „„ 18 „„„ iq 1989 2n 9585 2 , 32785323 22 
P(4, U )a 4 = - M 16 + 81?/ 18 - 297 U 19 + — u 20 -— M 21 + — M 22 

375123069 „ 2216152967 24 

ii H ii 24 , (58) 

10240 40960 ' v ; 

and for all N T > 5 up to 0(tt 8 ) in the correction 

P{N Tl u)a 4 = A u ^ t |^ t J 1 + 12 ^ tU 4 + 42 ^ tU 5_ 115|3 ^6 

597663 Ar 7 / 2 72206061 Ar \ 8 " 

7V T u 7 + 83 N T 2 H N T u 

2048 V 40960 j 

1 + 12 N T u 4 + 30 7V T u 5 - N T u® 

- 180 JV T u 7 + (^83 7V T 2 + ^ u 8 j . (59) 

Apart from the case N T = 1 all series are up to 0(u 8 ). For this case the number of contributing 
graphs was too large to give a reliable series expansion beyond 0(u 6 ). 
In the limit TV — > oo we find 

P(l,w)a 4 = 6u 4 + 44u 6 + 495w 8 + 5336w 10 , 
P(2,u)a 4 = 3 U 8 + 94 U 12 + 168 U 14 + ^ U 16 , 

At\A 



P(3,u)a 4 = 2u 12 + 72u w + ^u w + &3&u 2{ \ 
P{A,u)a 4 = V 6 + 72 U 20 + 120 U 22 + 1121 U 24 , 



-72 u 20 + 120u 22 + 1121u 24 , 
P(N T ,u)a A = ^ u 47V - (l + 12 A T u 4 + 20 7V T u 6 + (83 A r 2 - 147 7V T )u 8 ^. (60) 

Again we have less orders for the case N T = 1 . Interesting observations are the entirely positive signs 
and the vanishing of odd powers. 
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